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The phase transition of the Gross-Neveu model with A*' 
fermions is investigated by means of a non-perturbative evo- 
lution equation for the scale dependence of the effective av- 
erage action. The critical exponents and scaling amplitudes 
are calculated for various values of in d = 3. It is also 
explicitely verified that the Neveu- Yukawa model belongs to 
the same universality class as the Gross-Neveu model. 

The Gross-Neveu (GN) model [1] is one of the simplest 
models for interacting fermions. Nevertheless, in three dimen- 
sions our quantitative understanding beyond some universal 
characteristics of the phase transition has remined rather in- 
complete. The universality class of the GN-model in dimen- 
sions between 2 and 4 has been argued to be the same as the 
Neveu- Yukawa (NY) model [2] in 4 — e dimensions [3]. Both 
the large A'^ and e expansion indicate that a second order phase 
transition takes place for some critical value of the coupling 
constant if the number of fermion species A'' is larger than 
one [4]. The anomalous dimensions have been calculated up 
to the third order in the 1/N expansion [5], while some critical 
exponents have been computed to the order 1 /N in the phase 
with spontaneous symmetry breaking (SSB) [6]. In this letter 
we find the second order phase transition and calculate the 
critical exponents employing an analytical method based on 
nonperturbative fiow equations for scale dependent effective 
couplings. We directly obtain results for arbitrary dimension 
and without a restriction to large A''. Despite the presence 
of massless fermions we are able to investigate the symmetric 
phase. Due to the fermion fiuctuations the infrared physics is 
not trivial in the NY-language and requires a careful discus- 
sion of the critical exponents. Beyond the universal critical 
behaviour our method gives a description for arbitrary values 
of the GN-coupling away from the critical point. In particu- 
lar, we compute the non-universal critical amplitudes. 

The running couplings parameterize the effective average 
action Ft [7] which is a type of coarse grained free energy. 
It includes the effects of the quantum fluctuations with mo- 
menta larger than an infrared cutoff k. In the limit where the 
average scale k tends to zero Tt becomes therefore the usual 
effective action, i.e. the generating functional of IPI Green 
functions. In the limit ^ cxa it approaches the classical 
action. In a theory with bosons and fermions the scale de- 
pendence of Ffc can be described by an exact nonperturbative 
evolution equation [7,8] 

^r, [<^, V] = ^ Tr { (rf ) [0, ^] + nu)B^^^RuB 



dt 



RkF I 



(1) 



where t = ln(fe/A) with A some suitable high momentum 
scale. The trace represents a momentum integration as well 

(2) 

as a summation over internal indices and F), is the exact 
inverse propagator given by the matrix of second functional 
derivatives of the action with respect to bosonic and fermionic 
field variables. The infrared cutoff 




(27r)V(g-g') 



is parameterized by the bosonic and fermionic cutoff functions 
Rksiq) = q^Z„^krBiq), Rkpiq) = i /iZ^^trpiq)- We choose 



RkB 



Za.kq^ 



RkF — iZ^,k fl 



(2) 



where Za,k, Z^^k are wave function renormalizations. The 
momentum integration in Eq. (1) is both infrared and ul- 
traviolet finite. Equation (1) is an exact but complicated 
functional differential equation wich can only be solved ap- 
proximately by truncating the most general form of Fk ■ Once 
a suitable nonperturbative truncation is found the ffow equa- 
tion can be integrated from some short distance scale A, where 
Fa can be taken as the classical action, to fc — > thus solving 
the model approximately. 

The GN model is described in terms of a 0{N) symmetric 
action for a set of A'' massless Dirac fermions. The classical 
Euclidean action is given by 



Son = d X 



(3) 



(Here and in the following we distinguish with a bar the di- 
mensionful couplings.) The (pseudo)-scalar i7{x) is an aux- 
iliary non dynamical field which can be integrated out from 
the partition function, leading to the replacement (j{x) 
Gtp{x)tl}{x). Its vacuum expectation value ctq is proportional 
to the fermion condensate, ao = G < tptl^ >. The model 
is asymptotically free and perturbatively renormalizable in 
2 dimensions, hence it exhibits a non-trivial fixed point in 
d = 2 + e. It is 1/A'^ renormalizable in 2 < d < 4. 
The NY model whose classical action is 



Sny 



d X 



+ ^a\x) + y\x) 



(4) 
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has a Gaussian fixed point in d = 4 where it is perturbatively 
renormalizable and a non-trivial fixed point in d = 4 — e. Both 
models have, in even dimensions, a discrete chiral symmetry 
which prevents from the addition of a fermion mass term, 
while in odd dimensions a mass term is forbidden by space 
parity. Performing a large TV analysis the universal properties 
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of the two models arc argued to be the same in 2 < d < 4 
[3]:in such limit the two models are equivalent in the scaling 
region if we rescale ha to a and set G = jxr? . 

We consider a truncation of the effective action which 
contains a potential for the scalar field and a Yukawa term. 
In momentum space it reads {J dq = J d'*g/(27r)'') 



rk[a,ip,ij] = j (fxUk{rj) + j dq ^^a{-q)q^a{q) 
-Z^,ki)i{-q)i A'ip'iq) ~ j '^P hki)i{-q)(y{vW{<l - p) 



(5) 



The scalar potential is assumed to be a function of the invari- 
ant cr^(x) and we make the further simplification 



2 - 
TT t \ ''^k / 2/ X 2 N , flfc / 2/ N 2 n2 

Uk(<r) = ^(<T (a;) - aok) + -^{(t (x) - aok) 



1 ^k , 2/ X 



2 n3 



(6) 



The symmetric regime is characterized by the minimum being 
at ctqi. = 0. In the SSB regime a fe-dependent minimum af^f. / 
develops, whereas m| = 0. 

Inserting Eqs. (5), (6) into (1) we obtain a set of evolu- 
tion or renormalization group equations (RGE) for the effec- 
tive parameters of the theory in the two regimes. Integrating 
the RGE between some high momentum scale A and fc = 
we will find the phase transition point and extract the crit- 
ical exponents of the theory. We find it convenient to work 
with dimensionless quantities = Z^'^k'^~*h1 , gk = 
Z-^k^'^gk , bk = Z-^e-'-^k , Ck = Z-^k-^ml , p = 
^Zak^~'^cT^ , Kk = \Zak'^~'^alk, Uk = Ukk~'^ and we use 
^k = -gf etc. 

The evolution equation for the potential obtains from (1) 
by evaluating F^^' in the truncation (5) for a constant back- 
ground scalar field. We find 



k'i 



-VcrrB - ^yVB 



+ 2N' 



<+2K' + 2/(l + ^s) 



2hl~p + y{l + rpY 



Up) ■ (7) 



Here we have introduced the notation N' = 2^^^Af with 2"'^'^ 



the dimension of the 7 matrices and y = q /k 



we have defined v. 



2d+i^d/2p(.^^2). 



2^. Also 

dy 



In the SSB regime 



the evolution equations for the parameters Kk, gk and hk are 
then obtained as 



9tKfc = (2 -d-T]^)Kk 



[dpCk 



(8) 



+ -bkdtKk 
o 



(9) 
(10) 



gk 

dtgk = (2??CT + d - 4:)gk + 3 [d^Ck] ^ 
dtbk = (3r?a + 2d - 6)bk + 15 

L J Kk 

In the symmetric regime Kk = so we replace eq. (8) by 

dtCk = {r}a - 2)ek + [dpC,k]o ■ (H) 
The anomalous dimensions r)a and 77^ are defined as 

= —dt In Za,k , Vipik) = —dtlnZ^^k ■ (12) 



The wave function renormalizations Z parametrize the mo- 
mentum dependence of the propagators at zero momentum 
and (7 = aok ■ One finds 



?7a(fe) = 9a|^ J dyy"^^^ \{-^bkKk + gkf KkH(y,a^ 

- 2Nhl {2?yG{y, af - 2hlKkF{y, a)')] } _^ (13) 
(k) = 4/i^a„ I ^ j dy y'"^H(y, a)G{y, a) | (14) 
with 

H(y,a) = 2 FT N 7 vT 

Bk + ggkiik + y[{l + tb) - a(r)^rB + 2?/rs)J 

I + rp — ct{r)^rF + 2yrp) 



G{y,a) = 



F{y,a) 



y[l + rp - a{ri^rp + 2yrF)]'^ + 2/i|Kfc 
1 

y[l + rF - a{r]^rF + 2yfF)]^ + 2hlKk 



(15) 



Finally, the evolution equation for the Yukawa coupling ob- 
tains from taking derivatives of eq. (1) with respect to 4>,tp 
and a : 



dthl = (277V, +'na+d- 4:)hl 



4/ifc Vd dy y 
Jo 

2 (??,/, rF -I- 2yrF) 



d/2 



{narB+2yrB)H'\y,0)F{y,0) 



G\y,0)H{y,0) 
l + rp 



(16) 



The equations (13), (14) and (16) arc valid in both regimes 
provided we set Kk and appropriately. 

If we expand this set of equations in the coupling constants 
we recover the one-loop results obtained in the 4— e expansion 
for the NY model [2] 



dtg = -eg + (^g^ + 4Ngh' - 24Ar/i*) 



■yth' = -eh' + -^i2N + 3) 



4^ 



Vv = 



167r2 



(17) 
(18) 

(19) 



Moreover, after identifying the running coupling constant of 
the GN model as G — h\/ek we also recover the one- loop 
result obtained in the 2 + t expansion for the GN model [2] 



dtG = {d- 2)G -{N' -2)^+ 0{G^ 

271' 



(20) 



We numerically evolve the flow equations (8)-(14) from a 
large momentum scale A to fc ^ 0. The initial values of 
the parameters are chosen in such a way that Fa = Sgn- 
ZaA ^ 0, Z,pA = 1, hi = A, PA = 0, 5a = 0. Then 
sa = {Z^aGa)~^ is the only free parameter of the theory 
and plays the role of the temperature. This value has to be 
tuned in order to be near the second order phase transition. 
The value corresponding to the critical temperature Tc is de- 
noted by eAcr- For ZaA = 10"^", d = 3 and iV = 3 we find 
in our truncation eAcr = 1.878085212016 • 10® and the critical 
coupling in the GN model is GAcrA = 5.32. 
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The relevant parameter for the deviation form Tc is 5e — 
ga — EAcr = H(T — Tc) with constant H. In Fig. 1 we show 
the behaviour of the dimcnsionlcss couplings as functions of 
k for N = 3. As can be seen, they all reach a constant value 
in the symmetric regime (e^ > 0, = 0) corresponding to 
a nontrivial fixed point with vanishing beta functions. For 
(5e < the symmetry is broken, as expected: the coupling Ck 
goes to zero and the order parameter aok acquires a non zero 
value. 
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FIG. 1. Running couplings Sfc, -[g ^■'id Tmo' ^® func- 
tions of t = Infc/A given by dot-dashed, full, dashed and 
long-dashed lines respectively. 

In the following, whenever numerical results are reported, we 
fix d = 3 and N = 3. However, in Table 1. we summarize the 
results obtained for N = 2,4, 12. 

In the symmetric (high T) phase the fermions are massless. 
Their fluctuations induce a nontrivial dependence of and 
the renormalized scalar mass Tn%{k) — Z''^ {k)m^{k) on the 
scale k even away from the phase transition. This contrasts 
with the standard situation where the running of mR{k) es- 
sentially stops in the symmetric phase once k becomes much 
smaller than niR. The issue of critical exponents in a situation 
with two difi'erent infrared cutoffs k and mR is therefore more 
complex than usual. In a standard situation we would define 
the exponents 7 and z/ by following the temperature depen- 
dence of the unrenormalized and renormalized mass, m^(k) 
and m|j(fc) for fc ^ [9]. Here we define tfie renormalized 
mass at some fixed small ratio k/ruR by 



- 2 2 /, 1 2 
rriR = niRlkc) - turc 



.(fee), kc = rcTUR 



(21) 



with mjicrik) = efe^ on the critical trajectory. (In the numer- 
ical simulations we will fix the ratio Vc to be equal to 0.01). 
This mass corresponds to the only relevant parameter char- 
acterizing the critical behaviour. It is directly related to the 
deviation from the critical temperature T — Tc or 6e. In the 
following we use the arguments Scot friR interchangeably. We 
also define the inverse susceptibility or unrenormalized mass 
by 



= mjfZa{kc, TTifl). 



(22) 



Correspondingly, the critical exponents f and 7 are defined 
for fixed rn and we find 



1 91nm|j((5e) , 

1/ = - lim — - — — - = 1.041 

2 6e^o olnSe 

din fh^(5e) _ „„„ 
7 = hm ; ' = 1.323. 

Se-to alnde 



(23) 
(24) 



Ftom the definition (22) one has the relation 

din Z^{kc,rnR) 



21^ = j- 
A typical form of Za is 



din Se 



Zo 



rfiR + k^ 
A2" 



2 \ -2'''' 



k' 



and we conclude 

din Z„(kc., ffiR) 



d In 5e 



= -rjcr-u, -I = v{2-ria). 



(25) 



(26) 



(27) 



This is the usual index relation. We find fja — 0.729 and the 
index relation (27) yields 7 = 1.323, consistent with the value 
of 7 computed directly. The index 772, which vanishes in the 
standard situation, determines the dependence of ifiR on rc- 
For a more detailed understanding of the scale dependence 
we consider next the running of the renormalized mass and 
unrenormalized mass with k. We define 



f^(u g \ ^ 1 dln[m%{k,6e) - m|e^(fc)] 
1,^2 dt ''''' 



■y{k,5e) 



dln[m^{k, Se) 



r{k,Se)] I 



dt 



(28) 
(29) 



with m1^{k,5e) — m\cr{k.)Za{k,5e) . The relation m\ 
/Za implies the index relation 



^ = 20 -rja 



(30) 



which differs from the usual relation 7 = v{2 — r/^). For both 
k and ffiR sufficiently small and k » rfiR the indices i/, 7, rja 
and rj^ approach constant values independent of k and friR 
and we find 



£> = 0.502, 7 = 0.295, rj^ = 0.710, r]^ = 0.040. 



(31) 



These values agree well with Eq. (30) . In the opposite regime, 
k « rfiR, the running of the renormalized mass is only due to 
the anomalous dimension rja which is now different from the 
value (31). We find = 0.500, 7 = 0.000, r?a = 1.000, 7?^ = 
0.000. Again, these values agree well with Eq. (30) and the 
expectation 7 = 0. The nontrivial exponents O, rja in the NY- 
language correspond to the absence of rcnormalization effects 
for Gk for fc — > in the GN language. We note that for fixed 
Se the renormalized scalar mass (which corresponds to the 
inverse correlation length) scales as niR ~ k for tur « k and 
niR ^ \fk for niR, >> fc. The value of = 1 for k << tur, 
corresponds to r\2 in eq. (26). 

We next turn to the dependence on the deviation 5e at 
fixed k. We define the exponents as 



c X _ 1 ain[m|(fc,(5e) - mLr(fc)] , 
'^^'''^^^ " 2 d^e 



7(fc,(5e) 



dln\m}(k,&e) 



r{k, Se)] 1 



(32) 
(33) 



dlnSe 

In the limit Se ^ wc find i) = 0.520, 7 = 1.326. The 
indices (28), (29) and (32), (33) are related to the definition of 
the critical exponents (23), (24) by 



v = lim (^(fcc, Se) + ^(kc, Se)^) 
7 = lim (7(fcc, Se) + ^{kc, Se)^) 

<5e— »0 



(34) 
(35) 
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Here fcc(5e) is given by evaluating (21) at fixed Tc. Equations 
(34) and (35) yield respectively v = 1.041 and 7 — 1.326, 
consistent with the results of (23) and (24). 

We also have computed the (nonuniversal) critical am- 
plitudes which describe the dependence of fhR and m on 
the coupling Ga of the GN model. Observing Se/cAcr = 
GAcrS{l/GA) we obtain, for small deviations from criticality, 



ruR = Ai, 



SGa 



Gac 



- 2 
m 



SGa 



Gac 



(36) 



Our numerical values of the amplitudes A^, A^ are A^l K = 
0.016, A^jt^ = 0.212. 

We finally have investigated the low temperature phase 
with spontaneous symmetry breaking. Here the running of 
CTo stops for small k and the complications of the symmet- 
ric phase are absent. The critical exponent (3 is defined with 
(To = limfc^o cofc 



„ 1 ,. Ino-Q 
2 6e^o Inoe 



such that 



(To = Aa 



SGa 



Gac 



(37) 



(38) 



We find Afjjy/li = 0.008, fi = 0.903, in good agreement with 
the scaling relation (3 = ^{d — 2 + ija). In Fig. 2 we plot the 
condensate cro as a function of Ga- 







z 


Q 


A 
4 


1 9 
iZ 




u.yoi 


1 O/ll 


I.UIU 


i.UZO 


7 




1 '^9'^ 
l.OZo 


1 99S 


1 07f^ 
l.U ( 


J/^Z — rja ) 




1 ^9^ 


1 9'^0 
i.ZoU 


1 07f^ 
l.U ( 


n 
P 


74^ 
u. 1 ^0 




Q1 


QQR 




u. / ou 




QO'^ 


QQ1 


A IK 






OOQ 


01 4 


A / 


049 


919 


9'^'^ 
u.zoo 


QfiJ^ 
u.yuo 


A U/~\ 


n 007 
u.uu 1 


008 
U.UUo 


OOt^ 
U.UUO 


007 
U.UU / 


Vt 


0.561 


0.710 


0.789 


0.936 




0.066 


0.040 


0.027 


0.007 




0.541 


0.729 


0.765 


0.971 


GAcv 


9.989 


5.325 


3.613 


1.006 



TABLE 1. Critical exponents and amplitudes for different 
values of N 

The case N = 1 appears to be different from A'^ > 1 . We find 
a phase transition. For small Ga (large ca) eq. (20) is valid 
(N' = 2) and Gk scales according to its canonical dimension. 
The model is in the symmetric phase. For Ga > GAcr, GAcr = 
19.416 the mass term at the origin of the potential becomes 
negative, indicating spontaneous symmetry breaking. We find 
no scaling solution, neither for > nor for Kk > 0. This 
may suggest a first order transition. 

In conclusion, a simple truncation of the exact flow equa- 
tion for the effective average action gives a consistent picture 
for a second order phase transition for the GN model with 
Af > 2 in three dimensions. Wo have computed critical ex- 
ponents and amplitudes and we relate directly physical ob- 
servables like the correlation length or the order parameter 
to the value of the coupling G. By choosing different initial 
conditions we have also explicitoly verified that the Noveu- 
Yukawa model belongs to the same universality class as the 
Gross-Neveu model. The universal exponents are indepen- 
dent of the initial parameters, but not quantities like ao, the 
renormalized mass and the corresponding amplitudes. 
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FIG. 2. Formion-antifermion condensate (to/\/A as a 
function of the coupling Ga, in the range 

[3/4GAcr,3/2G Acr I- 

In the Table we report our results for different values of N. 
For all A'^ > 2 the existence of a second order phase transition 
is confirmed by our analysis. As can be checked, the scaling 
relations are well verified. To compare with existing results 
obtained in the 1/N expansion, let us fix AT = 12. In [6] the 
critical exponents have been calculated to the order 1/N: 



1 + 



3Ar7r2 



1.022, 7 = 1 -h 



8 



0- 



1 



ATtt^ 
16 
3Ar7r2 



1.068, 



: 0.955. 



(39) 



In the same paper Montecarlo simulations for A'^ > 12 are also 
reported. Conformal techniques have been used to calculate 
the anomalous dimensions to 0(1/A''^) [5] and yield, for A'' = 
12, r)^, = 0.013, rja = 0.913. However, such techniques rely on 
being exactly at the critical point and hence cannot be used 
to calculate the other exponents. 
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